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Abstract 

The analysis of new experimental data, obtained at the Proton Linear Accelerator 
of INR, with the aim to search for supernarrow dibaryons in the pd — * ppX\ and 
pd — > pdX2 reactions is presented. Narrow peaks with an experimental width of 5 
MeV at masses of 1904±2, 1926±2, and 1942±2 MeV have been observed in missing 
mass MpXi spectra. In the missing mass Mx 1 spectra, the peaks at Mj x = 966 ± 2, 
986±2, and 1003±2 MeV have been found. The analysis of the data obtained leads to 
the conclusion that the observed peaks in M p x 1 spectra are most likely supernarrow 
dibaryons, the decay of which into two nucleons is forbidden by the Pauli exclusion 
principle. An alternative interpretation of the spectra by assuming a decay of the 
supernarrow dibaryons in "exotic baryon states" with masses Mx 1 is discussed. 
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1 Introduction 



The experimental search for dibaryons continues since over 20 years (for recent reviews see 
The dibaryons considered couple usually to two nucleons and reflect their symmetry. 
Therefore, such dibaryons have estimated decay widths from several MeV up to a few 
hundreds MeV. However, their production cross sections in reactions involving nucleons 
and pions are small and the spectra are dominated by a large background. This has lead 
to contradictory results and the dibaryons decaying directly into two nucleons are not 
unequivocally established. 

We consider a new class of dibaryons - supernarrow dibaryons (SNDs), the decay of 
which into two nucleons is forbidden by the Pauli exclusion principle . Such dibaryons 
with masses M < 2m^ + m n , where mjv(m 7r ) is the nucleon (pion) mass, can decay into 
two nucleons, by emitting a third particle, normally a photon. The decay widths of these 
dibaryons are predicted to be < 1 keV ||. 

In the framework of the MIT bag model, Mulders et al. B calculated the masses of dif- 
ferent dibaryons, in particular, iViV-decoupled dibaryons. These dibaryons D(T = 0; J p = 
0~, l-,2-; M = 2.11 GeV) and D(l; 1~; 2.2 GeV) correspond to the states 13 Pj and 31 P X 
forbidden in the NN channel. However, the masses of these dibaryon exceed the pion pro- 
duction threshold. Therefore, these dibaryons can decay into the -kNN channel and their 
decay widths are larger than 1 MeV. 

Using the chiral soliton model, Kopeliovich |7|] predicted that the masses of D(T = 
1, J p = 1 + ) and _D(0,2 + ) iViV-decoupled dibaryons exceeded the two nucleon mass by 60 
and 90 MeV, respectively. These values are lower than the pion production threshold. 

In the framework of the canonically quantized biskyrmion model Krupovnickas al. |5| 
obtained an indication of the existence of one dibaryon with J=T=0 and two dibaryons 
with J=T=1 with masses smaller than 2mjy + rn n . 

All these results are model calculations of very limited predictive power and provide just 
a motivation to look for SNDs. In the following we summerize the experimental attempts 
to look for such dibaryons so far. 

In ref. H] dibaryons with exotic quantum numbers were searched for in the process pp — > 
pp'j'j. The experiment was performed with a proton beam from the JINR phasotron at an 
energy of about 216 MeV. The energy spectrum of the photons emitted at 90° was measured 



and showed two peaks. This behavior of the photon energy spectrum was interpreted as 
a signature of an exotic dibaryon resonance with a mass of about 1956 MeV and possible 
isospin T = 2. However, the results of ref. jllj make the possibility of a production of 



dibaryons with T = 2 in this reaction questionable. So, additional careful studies of the 
reaction pp — > pp^ are needed to more correctly understand the nature of the observed 
state. 

An analysis of the Uppsala proton-proton bremsstrahlung data [[Llj] looking for the 
presence of a dibaryon in the mass range from 1900 to 1960 MeV gave only upper limits 
of 10 and 3 nb for the dibaryon production cross section at proton beam energies of 200 
and 310 MeV, respectively. This result agrees with the expected values of the cross section 
obtained at the conditions of this experiment in the framework of the model of dibaryon 



production and decay suggested in ref. |15| and does not contradict the data of ref. 

In our previous studies [j^ -|5] we investigated the reaction pd — ► pX in order to search 
for SNDs. In this process SNDs can be produced only if the nucleons in the deuteron 
overlap sufficiently, such that a six-quark state with deuteron quantum numbers can be 
formed. Then, an interaction of a meson or another particle with this state can change its 
quantum numbers so that a metastable state is formed. The experiment was carried out at 
the proton beam of the Linear Accelerator of INR using the two- arm spectrometer TAMS. 



As was shown in ref. |JT4||lo|l, the nucleons and the deuteron from the decay of SND into 
'-fNN and jd have to be emitted into a narrow angular cone with respect to the direction 
of the dibaryon. If a dibaryon decays mainly into two nucleons, then the expected angular 
cone of the emitted nucleons is about 50°. Therefore, a detection of the scattered proton 
in coincidence with the proton (or the deuteron) from the decay of particle X at correlated 
angles allowed to suppress effectively the contribution of the background processes and to 
increase the relative contribution of a possible SND production. As a result, two narrow 
peaks in missing mass spectra have been observed at M =1905 and 1924 MeV with widths 
equal to the experimental resolution (~ 3MeV) and with 4.8 and 4.9 standard deviations 
(SD), respectively. The analysis of the angular distributions of the charged particles (p 
or d) from the decay of particle X showed that the peak found at 1905 MeV most likely 



corresponds to a SND with isotopic spin equal to 1. In ref. [|15[] arguments were presented 
that the resonance at M =1924 MeV could be SND, too. 



2 Experiment 



In this paper we present a new study of the pd — > pX reaction at the Linear Accelerator 
of INR with 305 MeV proton beam using the spectrometer TAMS where the proton and 
deuteron from the decays X — > pX 1 and X — > <iX 2 are measured in coincidence with the 
scattered proton. 

In the experiment, CD 2 and 12 C were used as targets. The scattered proton was detected 
in the left arm of the spectrometer TAMS at the angle 9 L = 70°. The second charged particle 
(either p or d) was detected in the right arm by three telescopes located at Or = 34°, 36°, 
and 38°. 

A trigger was generated by four-fold coincidences of the two AE detector signals of 
the left arm combined with those of any telescope of the right arm. The events contained 
information about time-of-flights and full energies of two particles detected in coincidence 
in the left and right arms of the spectrometer. The energy resolution was 4 MeV (7 MeV) 
for the left (right) arm. A time-of-flight resolution better than 0.5 ns was achieved. Each 
valid event including two time-of-flights and two energies were stored event by event and 
then analyzed off line. 

An off-line identification of different particles was performed by means of their energies 
E and time-of-flights t. In this way, protons and other charged particles were identified via 
characteristic loci observed in two dimensional diagrams of time-of-flight versus energy. As 
an example such an experimental E — t distribution for the right arm detector at Or = 38° is 
displayed in Fig. 1. This figure shows the band of proton events, presenting the dependence 
of the proton time-of-flight on its energy, and a spot corresponding to the deuteron from 
the elastic pd scattering. 

Several software cuts have been applied to the mass spectra. The energy of the scattered 
proton was limited by an interval of 50 < E 2 < 150 MeV, that corresponded to the interval 
of measured dibaryon masses 1980 > M pXl (M dX2 ) > I860 MeV/c 2 . In order to suppress 
the background from the reaction p 12 C — > dY, where Y stands for the unobserved particles, 
in Mdx 2 spectra, we omitted events with a deuteron energy higher than the energy of the 
elastically scattered deuteron in the pd — > pd reaction at the given angles. For the M pXl 
spectra the energy interval for the proton from the SND decay was determined by the 
kinematics of the decay of X into the 7iVW channel. This interval was equal to 50 < 



Ep < 100 MeV. Such a cut allows to suppress essentially the contribution from the reaction 
pd — > ppn. The main contribution to this reaction at the energy under consideration is given 
by one-nucleon pole diagrams when one of the nucleons is a spectator. These diagrams give 
the maximum contribution when the kinetic energy of the spectator nucleon approaches 
zero. 

Missing mass spectra of M v x x and Max 2 were determined with the expression 

Kx^dx,) =m 2 d + 2ml + 2m ^ E \ ~ E f) ~ 2E \ E ) + 2 ™f cos 6 L (1) 

with the additional condition that the proton from the X — > pX\ and the deuteron from 
the X — > dX 2 decay were detected in the right arm detector of the TAMS spectrometer. In 
this expression m p (d) is the proton (deuteron) mass, E\ and Ej are the total energies and 
Pi and pf are the momenta of the incident and scattered protons, respectively. 

The missing mass M d x 2 spectrum of the reaction pd — > pdX 2 for the angle 9r = 38° (Fig. 
2) has a peak at the deuteron mass. This peak corresponds to the elastic pd scattering. 
A measurement of this reaction at different angles of the right and left arms was used to 
calibrate the spectrometer. The over all mass resolution of the spectrometer was ~ 5 MeV, 
the angular resolution ~ 1°. 



3 Results and Discussions 

We first consider the pd — > pX — > ppXi reaction. Figs. 3a-3c depict the experimental 
missing mass M p x 1 spectra obtained with the CD 2 and 12 C targets, where (3a), (3b), and 
(3c) correspond to a detection of the second proton in the right arm detector at 9r = 34°, 
36°, and 38°, respectively. The background in these spectra are interpolated by polynomials. 
Three peaks at M pXl = 1904 ± 2, 1926±2, and 1942±2 MeV are clearly visible in these 
spectra. The first two of them confirm the values of the dibaryon mass obtained by us 
earlier [[12] - 15 1 , whereas the resonance at 1942 MeV is a new one. This additional resonance 



is due to the different kinematical coverage of the new experiment. In this experiment we 
could decrease the detection threshold of the scattered proton down to 50 MeV giving us a 
broader acceptance in the missing mass. 

The experimental missing mass spectra, obtained with the carbon target, are rather 
smooth. This smoothness is caused by both, an essential increase of the contribution 



of background reactions in the interaction of the proton with the carbon and the Fermi 
motion of the nucleons in the nucleus. The latter increases essentially the angular cone of 
the emitted nucleons. In consequence, it is not possible to see peaks of SNDs on the carbon 
target. As the experiment with the carbon target resulted in the rather smooth spectra, 
all structures, appearing in the experiment with the CD 2 target, have to be attributed to 
an interaction of the proton with the deuteron. 

The experimental spectra in Figs. 3a-3c are compared with the prediction of the the- 
oretical model of SNDs D(T = 1, J p = 1^) production constructed in the one meson 
exchange approach |0| and normalized to the values of peaks in Fig. 3a. Our calculation 
for the isovector SND D(T = 1, J p = 1 ± ) with the mass M=1904 MeV showed that the 
contributions of such a dibaryon to the spectra at angles 34°, 36°, and 38° must relate as 
1 : 0.92 : 0.42. For the isoscalar SND J D(0,0 ± ) one obtains 1 : 0.95 : 0.67. The biggest 
contribution of the SND with M = 1926 MeV is expected at 6 R = 30° H- 34°. However, 
the angles 30° -r- 33° were not investigated in this work. The calculation of the ratio of the 
contributions to the missing mass spectra at the angles 34°, 36°, and 38° gave 1 : 0.85 : 0.34 
for T = 1 and 1 : 0.71 : 0.46 for T = 0. 

Nucleons from the decay of the SND with M = 1942 MeV have a wider angular distri- 
bution due to the the higher mass and consequently the higher transverse momenta in the 
decay X — > pX\ according to the calculations with a maximum in the region of 26° -j- 32°. 
In the region of the angles under consideration here, the contributions of the D(l, l 1 * 1 ) are 
expected to behave as 1 : 0.6 : 0.2. For the SND D(0, ± ) we have 1 : 0.78 : 0.55. 

All these predictions for the SNDs are in agreement with our experimental data within 
the errors. However, the analysis of the reaction pd — > ppX\ only in the considered angular 
range does not allow to determine an isotopic spin of the SNDs. 

If the observed states are iViV-coupled dibaryons decaying mainly into two nucleons 
then the expected angular cone size of emitted nucleons must be more than 50°. Therefore, 
their contributions to the missing mass spectra in Fig. 3a-3c would be nearly constant and 
would not exceed a few events, even assuming that the dibaryon production cross section 
is equal to that of elastic pd scattering (~ 40^xb/sr). Hence, the peaks found most likely 
correspond to SNDs. 

The missing mass Mdx 2 spectrum of the reaction pd —> pX — > pdX 2 , for the sum of 
angles 9r = 34° and 36° is shown in Fig. 4. As seen from this figure, the reaction pd — > pdX 2 



gives a very small contribution to the production of the dibaryons under study. 

On the other hand, it is expected n,|T3] that isoscalar SNDs contribute mainly in the 
'yd and isovector SNDs in the -yNN channels. As the main contribution of the indicated 
dibaryons is observed in pX x channel, it is possible to suppose that X x — > and all found 
states are isovector SNDs. A more precise conclusion about the value of the isotopic spin 
of the observed SNDs could be obtained by the study of the reaction pd —> npX\. The 
experiment on the SND photoproduction in the process 7<i — ► ti ± D with polarized photons 
could separate the D(T = 1, J p = 1 + , S = 1) and D(l, 1~, 0) states [[Uj where S is the spin 
of the SND. 

The sum of missing mass spectra of the reaction pd — > ppXi over angles 9r = 34° 
and 36°, where the contribution of the SNDs is maximum, is presented in Fig. 5a. This 
spectrum was fitted by a polynomial for the background plus Gaussians for the peaks. The 
number of standard deviations (SD) is then determined from this spectrum as 

"7^2= (2) 

where N e ff is the number of events above the background curve and A^acft is the number of 
events below this curve. Taking nine points for each peak, we have 6.0, 7.0, and 6.3 SD for 
the resonances at 1904, 1926, and 1942 MeV, respectively. The widths of these resonances 
are equal to the experimental resolution of ~ 5 MeV. 

An additional information about the nature of the observed states has been obtained 
by studying the missing mass Mx x spectra of the reaction pd — > ppX\. If the state found 
is a dibaryon decaying mainly into two nucleons then X\ is a neutron and the mass Mx x is 
equal to the neutron mass m n . If the value of Mx 1 , obtained from the experiment, differs 
essentially from m n then Xi — > 772, and we have the additional indication that the observed 
dibaryon is SND. 

The simulation of the missing Mx x mass spectra for the reaction pd — > ppX\ has been 
performed assuming that the SND decayed as SND — > jd( 31 So) — > "ypn through the two- 
nucleon singlet state 31 Sq ||[15|]. Such a decay is characterized in the rest frame by a narrow 
peak near the maximum photon energy in the probability distribution of the dibaryon decay 
over an emitted photon energy. As a result of this simulation, in the missing Mx x spectra 
three narrow peaks at Mx t = 965, 987, and 1003 MeV have been predicted. These peaks 
correspond to the isovector SNDs with the masses 1904, 1926, and 1942 MeV, respectively. 



The result of this simulation without the influence of the detector resolution is shown in 
Fig. 6. Fig. 6a shows the Mx 1 mass spectrum when the proton from the SND decay 
SND — > -ypn is emitted into the whole allowed angular region. The acceptance of the 
detectors cuts into the widths of these peaks. The missing Mx x mass spectra, taking into 
account the acceptance of the detectors at the angles 34° and 36°, is presented in Fig. 6b. 
The experimental resolutions widen these peaks to 6 — 8 MeV. 

Fig. 5b depicts the sum of the missing mass Mx x spectra for the angles 6% = 34° and 
36°. As is seen from this figure, besides the peak at the neutron mass, which is caused by 
the process pd —> ppn, peaks are observed at 966 ±2, 986 ±2, and 1003 ±2 MeV. Hence, for 
all states under study, we have Xi — > 'jn in support of the assumption that the dibaryons 
found are SNDs. 

It should be noted that the peak at Mx x = 1003 ± 2 MeV corresponds to the resonance 



found in ref. Jl7j and was attributed to an exotic baryon state N* below the nN threshold. 
In that work, the authors found altogether three such states with masses 1004, 1044, and 
1094 MeV. This is not in contradiction to our interpretation, since in principle SND could 
decay according to X — > pN*. The possibility of the production of iViV*-coupled dibaryons 
was considered in [J2J. 

If the exotic baryon states decay into N* — > then they would contribute to the 
Compton scattering on the nucleon. However, the analysis [|18j of the existing experimental 
data on this process excludes N* as intermediate states in the Compton scattering on the 
nucleon. 

In ref. |H| it was assumed that these states belong to totally antisymmetric 20-plet of 
the spin-flavor SU(6)fs group. Such a N* can transit into a nucleon only if two quarks 
from the N* participate in the interaction ||20|| . Then the simplest decay of iV* with the 
masses 1004 and 1044 MeV would be N* — > 77^. This conjecture could be checked, in 
particular, by studying the reactions 'jp —>■ 'jX or 7p — > nX at a photon energy close to 
800 MeV. 

Taking into account the mentioned connection between the SNDs and the resonancelike 



states X\, it is possible to assume that the peaks, observed in [[IT]] at 1004 and 1044 MeV, 
are not exotic baryons, but they are the resonancelike states X v — > caused by the 
existence and decay of the SNDs with the masses 1942 and 1982 MeV, respectively. Such 
X v would be not real resonances and cannot give contribution to the Compton scattering 



on the nucleoli. 



4 Conclusions 

The result of the study of the reaction pd — > ppX\ is that three narrow peaks at 1904, 
1926, and 1942 MeV have been observed in the missing mass M pXl spectra. The analysis 
of the angular distributions of the protons from the decay of the pXi states and the data 
of the reaction pd — > pX — > pdX 2 showed that the peaks found can be explained as a 
manifestation of the isovector SNDs, the decay of which into two nucleons is forbidden by 
the Pauli exclusion principle. The observation of the peaks in the missing mass M Xl spectra 
at 966, 985, and 1003 MeV is an additional confirmation that the dibaryons found are the 
SNDs. 
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5 Figure captions 

1. Energy E against time-of-flight t scatter plots for the right arm detector at Or = 38°. 

2. The missing mass M dX2 spectrum for Or = 38°. 

3. The missing mass M pXl spectra obtained with CD 2 (the points with the statistical 
errors) and 12 C (the bars) targets; (a) R = 34°, (b) R = 36°, (c) R = 38°. The 
solid curves are normalized theoretical predictions. The dashed curves correspond to 
the background interpolated by polynomials. 

4. The sum of missing mass M dX2 spectra for the reaction pd — > pX — > pdX 2 for the 
angles 9 R = 34° and 36°. 

5. The sum of missing mass M pXl (a) and M Xl (b) spectra for the angles Or = 34° and 
Or = 36°. The dashed and solid curves are results of a fit by polynomials for the 
background and Gaussians for the peaks, respectively. 

6. The result of the simulation of the missing M Xl mass spectra; (a) - the proton from 
the decay X — > pX\ is emitted in the whole allowed angular region, (b) - this proton 
is detected at Or = 34° and Or = 36°. 
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